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The missing of a Keplerian fall-off in the observed galaxy rotation curves represents classical 
evidence for the existence of dark matter on galactic scales. There has been some recent activity 
concerning the potential of modelling galactic systems with the help of general relativity. This 
was motivated by claims that by the use of full general relativity dark matter could be made super- 
fluous. Here we focus on possible axisymmetric and stationary solutions of Einstein's equations 
with rotating dust. After a short review of the current debate we pursue the idea of approaching 
such relativistic models in a Newtonian language. We analyse rigidly as well as differentially 
rotating Newtonian and Post-Newtonian spacetimes and find that it is necessary to incorporate a 
Post-Newtonian term in order to make physical sense. 
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1. Motivation 

The evidence for dark matter stems from various observations at various scales but it is not 
necessary that the effects have all the same origin. This paper is motivated by the galactic dark 
matter problem: in contrast' to the expectation from Newtonian dynamics the tangential velocities 
of spirals are found to be approximately flat and not decreasing up to large radii. 

The standard resolution lies in postulating a large halo of dark matter around galaxies. In 
contrast to this Cooperstock and Tieu (CT) recently suggested that the proper modelling of a galaxy 
through full general relativity (GR) explains the findings without dark matter [Q]. Their galaxy 
model is a stationary axisymmetric (comoving dust) solution based upon the following metric 

As^=e''-'''{d? + dz^)+r^e-"d(p^-e''{dt+Nd(pf . (1.1) 

Even though fields and velocities in the galaxy are small, according to CT, the Newtonian approxi- 
mation would not be valid for a galactic system that is under the influence of gravity alone because 
of intrinsic non-linear contributions of non-negligible size from GR. The GR model would then 
be able to reproduce realistic density profiles after fitting flat rotation curves. The according total 
galactic mass comes out of modest value, c.f. [^, that is smaller than the newtonianly required 
value and larger than values from Modified Newtonian Dynamics (MOND). 

By now various authors have pointed to problems with the consistency of the CT model. In 
the Komar mass within the CT model is analysed and it is found that the energy momentum tensor 
becomes ill-behaved at z = 0, revealing the existence of an additional exotic matter source. Thus, in 
the CT model, the location of exotic matter would only be transfered to the galactic plane, a finding 
also supported in [Q]. In it is argued that Post-Newtonian corrections should already enfold 
non-linear effects if they are present and that corrections should intrinsically be small in the weak 
field limit. In ||^ an inconsistency between the CT comoving frame and the condition of differential 
rotation is revealed; if the inconsistency is removed the rotation curve became ordinarily Keplerian. 
It is argued in that the CT model gives wrong predictions for the lateral density profile of the 
Milky Way and for its local density. On the other hand, in it is found that the CT model seems 
to hide infinite mass at large distances. Generally it is shown in |§] that interior stationary and 
axisymmetric solutions - like CT - cannot be continued into according exterior solutions in order 
to become global solutions without evoking singularities. 

It is interesting to note that viable disk models for galactic systems with differential rotation do 
not exist in GR by now. For the case of rigid rotation [|l^] an exact solution is known which has been 
obtained via the inverse scattering method; already its formulation is technically rather elaborate. 
Though problematic, the CT hypothesis has initiated further alternative ideas on galactic rotation 
curves. Using a similar GR model in [O] it is found that the amount of dark matter needed is 



reduced by 30%. In [ |12| ] an axisymmetric disk solution of the Einstein equations in six dimensions 
is constructed. The modification of gravity via extra dimensions then leads to flat rotation curves if 
orbits within the disk are stable. 

Here we want to compare the CT model with a rotating (Post-) Newtonian model. For this we 
first have to ask some general questions on axially symmetric and stationary solutions in GR. 



'in Newtonian gravity a flat rotation curve is possible, see Mestel's disk But the price is an infinite total mass. 
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2. Rotating (Post-) Newtonian Spacetimes 



In notation (x°,x\x^,x^) = {t, r, (p,z) and signature (— , +, +, +) the most general axisymmet- 
ric {dipg^y = 0) and stationary {dtg^v = 0) spacetime in GR with the four free metric functions 



U {r,z),k{r,z),W {r,z),A{r,z) is named after Lewis and Papapetrou (LP) and reads [13| 



d/ 



e^''{dr' + dz.')+W'd(p^ - e'"" {dt + Ad(p) 



(2.1) 



We can simplify ( [2. ID a bit more, but only under crucial assumptions. If and only if the metric 
function W is harmonic it can be transformed toW = r. Consider a complex coordinate transfor- 
mation f{r + iz) = W + iV introducing an additional potential V . Then we have from p = W{r,z) 
and h = V{r,z) the differentials dp = ^dr + ^dz and dh = ^dr + ^dz. After inserting into 



(2.1), written in terms of p,h, and requiring formal invariance as compared to the original metric 
we see that the mixing terms have to vanish. That is exactly provided by the Cauchy-Riemann 
equations for W and V: 



dw 
IF 



f and ^ 



^ . Moreover, with the help of these equations, we 



see that the coefficients of dr and dz can be combined and so we reobtain ( |2.1| ). 

Thus we have shown that it is possible to simplify the general LP form (11) by allowing 
W = r, which is only possible if W (and V) be harmonic^: A^^^^VK = 0. Upon this constraint we can 
write down the LP metric in isotropic coordinates (or Weyl gauge) 



ds' 



-2U 



e^\dr^ + dz^) + r^d(p^ - {dt +Ad(pf . 



(2.2) 



The spacetime applied in the CT model was ( |1.1[ ). Obviously, the CT metric does not belong to the 
class of the most general stationary and axisymmetric spacetimes; it only belongs to the subclass 
of LP solutions in the Weyl gauge, and is therefore less general. 

Now we ask what solution could potentially be a Newtonian counterpart to the CT model ([Ell). 
The 'Newtonian approximation' [[T^] that is the metric that reproduces Newtonian physics is 



d^" 



- ( 1 + 2(/) )df2 + dr^ + r^dcp^ + dz^ , 



(2.3) 



where ^ = ${r,z) is the Newtonian gravitational potential. For simplicity, we start with rigid 
rotation in ( p3| ) via (p = (p' — cot and CO =const. The exact result can be brought to the form 

-I 2 



ds^ = {dr^+dz^) + - 



1+20 



CO^r^ 



-r^dtp^ - 



(1+20- 



dt + 



(l+2(/)-w2r2) 



d(p 



(2.4) 



(1+20 

In this form we can directly compare it with the LP metric in Weyl gauge (2.2), and we notice a 
discrepancy at linear order in , looking at the d(p^ term. Interestingly, the rigidly rotated Newton 



metric ( |2.4| ) is not consistent with the Weyl subclass of the LP solution ( |2.2[ ). 

One could now ask whether the situation might be easily cured with the help of a coordinate 
transformation. We show that this is not possible. Above we have derived the exact conditions 
under which the general and the isotropic LP metric can be transformed into each other, that is the 
function W must fulfil the two-dimensional Laplace equation A^^)^ _ q jj^ (j^g ^^^g of (2.4) we 
identify W = + 20. Expanding to linear order and applying the Laplacian yields 



A(2)w = rA(^'0 + 0,,- = 47rGpr + 0,,- 



(2.5) 



^This condition for W holds for exterior (vacuum) solutions that are stationary and axisymmetric, c.f. [p^p^] 
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after using the Poisson equation. We will now show that A^^VK does not vanish in general. 

Given the problem of solving the Laplace equation with boundary conditions for a disk-like 
distribution of matter, the solution for the potential can be obtained [ [l7| ] via separation of variables: 
(^(r,z) = S{k)jQ{kr)e^^^^^dk. A surface mass density r(r) is then characterised by its Hankel 



transform S{k) = —2kG Jq Jo{kr)'L{r)rdr . We use these expressions for the evaluation of (2.5). 

Case (I) z / 0. Outside the disk the Newtonian potential fulfils the Laplace equation, such that 
( 2^ ) becomes A^^^ W = — /q°° S{k)J\ {kr)ke^'^^~^ dk which will not vanish in general. We demonstrate 



this e.g. with the Mestel model which is characterised by a surface mass density that falls off 
inversely with the distance L{r) = Loro/r and Hankel-transforms as S{k) = —2nGLoro/k. Using 
this we can integrate and obtain A^-^^W = iKGLoro (^j — —0=^ at z / 0. 

Case (II) z = . We want to show that ( ^ ) is non-zero also here. Let us assume the contrary 
and see what happens. If we assume that A^^^VK = was true then equation (2.5) gives an identity. 



This we integrate over z for some £ > and then revoke the operation with the appropriate limit, 
—AnGrlimg^o /^g 5(z)r(?")dz = limg^o J-e lo S{k)Ji {kr)ke^'^\"\ dkdz . Since the exponential term 
on the right hand side serves as a damping factor, the modulus of the integrand will reach its 
maximum at z = 0. Thus, as an upper estimate, we can set the integrand of the right hand side to be 
constant in z and therefore the integration and limit procedure give zero. Then, for all other z the 
expression will be zero more than ever and we have 47iGL{r)r = 0: this will not hold generally for 
any realistic model, and thus A(^)W(r,z) = is not true at the surface z = either. 



Notably, the classical Newton metric (2.3) cannot be made compatible with the LP metric in 



Weylform and is thus also not compatible with the CT model. Additionally, as all stationary and 



axisymmetric vacuum solutions belong to the Weyl class [16], it follows that such a galaxy does not 
have a rotating classically Newtonian exterior solution. We go one step further and consider the 
Post-Newtonian^ (PN) metric with an additional PN potential = ^{r,z) , 

ds^ = -(l + 20)df2 + (l-2vA)(dr2 + r2d<p2 + dz^). (2.6) 

To include one more complication we consider the differential rotation of this metric via ft) = 
(o{r,z) and arrive at a lengthy expression. Unfortunately, the resulting metric exhibits direct time 
dependence in some coefficients and is therefore only reasonable in a strictly local sense. To 
preserve stationarity we approximate it by allowing for only small time intervals or equivalently 
for small angles of rotation. We find that for the simplest PN case, i.e. ^ = y, the differentially 



(and the rigidly) rotated metric is consistent - neglecting terms 0{(}>^) - with the Weyl form ( p!2| ). 

3. Summary and Conclusion 

Up to day, there does not exist an intuitive and applicable GR solution which could help to 



realistically model a galactic system - except for the complicated solution [10] with stiff rotation 



CT have proposed a solution but as we have seen arguments from various directions have revealed 



^Sometimes, this metric with = is referred to as the 'Newton metric' in the literature. The reason for this 
nomenclature might be that the order of magnitude of the coefficient of the spatial part and the order of the Newtonian 
correction are the same. However, conceptually this makes a considerable difference. In classical Newton Gravity there 



exists no curvature of space, the three-space is always euclidian. This is exactly reflected in the Newton metric (2.3) and 
therefore we refer to ( 2.6 ) as the PN approximation; for an extensive discussion see e.g. chp. 39 of [|l4]]. 
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unphysical features in the model. We have shown that the CT solution does not belong to the class 
of the most general axially symmetric and stationary solutions, the LP class. Therefore the CT 
solution is less general and this restriction might be connected with the problems of the model. 

A more intuitive understanding of the CT model and its breakdown would be instructive. 
For this we approached the CT model by Newtonian means. We found the interesting result that 
the classical Newton metric after rotation does not lie in the same class as the CT metric but is 
more general and cannot be simplified. Only the incorporation of a PN potential gives rise - after 
differential or rigid rotation - to equivalent metrics. Thus the proper weak field limit of such a 
galaxy must be Post-Newtonian. A complete analysis of the resulting dynamics via the Amowitt- 
Deser-Misner formalism is the concluding step and is going to be published elsewhere. 
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